Ultraslow light in inhomogeneously broadened media 
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We calculate the characteristics of ultraslow light in an inhomogeneously broadened medium. We 
present analytical and numerical results for the group delay as a function of power of the propagating 
pulse. We apply these results to explain the recently reported saturation behavior [Baldit et ai, 
Phys. Rev. Lett. 95, 143601 (2005)] of ultraslow light in rare earth ion doped crystal. 
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The usage of a coherent field to control the opti- 
cal properties of a medium has led to many remark- 
able results such as enhanced nonlinear optical effects 
1 , 01 , electromagnetically induced transparency (EIT) 
3|, lasine without inversion 0, 0, @i ultraslow light 
st orag e and retrieval of optical pulses [T^l 
and many others IT3 . ITEl . Most of these effects 
rely on quantum interferences which are created by the 
application of a coherent field. The coherent field opens 
up a new channel for the process under consideration. 
This interference effect produces the EIT dip or a hole in 
the absorption profile. The ultraslow light emerges as the 
EIT dip could be very narrow. It has been realized that 
in principle one could also use two level nonlinearities in 
presence of a strong pump. For a homogeneously broad- 
ened medium a hole can emerge if the transverse and 
longitudinal relaxation times are quite different. Under 
these conditions the hole has a width of the order of Ti 
and this is being referred to as the effect of coherent pop- 
ulation oscillation |l7j| . Bigelow et al. did experiments 
in this regime using ruby as the material medium which 
can be modelled as a homogeneously broadened system 
[TsL [T9| . Some studies on slowlight in inhomogeneous 
broadened medium exist pol [21i] . In an earlier paper 
the present authors had considered the case of inhomo- 
geneously broadened gaseous medium where the Doppler 
effect is important 21]. We considered the case of sat- 
uration absorption spectroscopy. This leads to the well 
known hole in the Doppler profile. The width of this hole 
was of the order of 1/Ti which is about two times 1/T2. 
In the inhomogeneously broadened gaseous medium the 
group index of the order of 10^ was obtained. The re- 
cent experiment of Baldit et al. reports group delays of 
the order of 1.1 s in rare earth ion doped crystal which 
has strong inhomogeneous broadening In this case 
all the relaxation times are quite different - Ti — 8 ms; 
T2 = 3 /is; inhomogeneous line width Tinh = 1-3 GHz. 
The width of the whole is essentially determined by Ti 
and hence one gets very large delays. Baldit et al. did 
present a theoretical model based on homogeneous broad- 
ening of the medium where as to obtain agreement with 
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FIG. 1: A schematic diagram of two level atomic system with 
ground state \g) and excited state |ei); The pump (cuc) and 
probe (oip) fields are co-propagating. 



experiments inhomogeneous broadening is to be included 
as alluded by them |23j . 

In this paper, we consider a system of inhomoge- 
neously broadened two level atoms interacting with co- 
propagating pump and probe fields. We use the well 
known susceptibility |24j and average it over the inho- 
mogeneous distribution to calculate the group index. We 
derive a number of analytical results and show how these 
can be used to understand the experimental results of 
Baldit et al. For example we show that in the limit of very 
small detuning of the probe from the pump the group de- 
lay goes as y/S for large S. The group delay also peaks 
at about S=0.9. The value of group delay increases as 
the detuning 5 increases. We further present detailed 
numerical results. 

In order to understand the experimental results of 
Baldit et al, we consider a two level system as shown 
in figure n Here we define all fields as 



E,iz,t) = £,{z,t) e-*(-'*-'=-) 



-t- c.c. 



(1) 



where £i is the slowly varying envelope of the field. The 
pump field at frequency iVc and the probe field at fre- 
quency ujp are co-propagating through the medium. The 
linear susceptibility x('-^p) is obtained by solving the den- 
sity matrix equations for the two level system of figure ^ 
that is by calculating the density matrix element p^g to 
the first order in the probe field but to all orders in the 
co-propagating pump field. The dynamics of population 
and polarization of the atoms in the two-level configura- 
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tion are given by 
^ 1 
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Pe 



-P., = 1 (2) 



where Ti and r2 are the longitudinal and transverse re- 
laxation times respectively. The density-matrix elements 
in the original frame are given by p^ge""^"*, p^^ , and p^^. 
The detunings A, S and the Rabi frequencies are defined 
by 

Iff f 7ff f 

= UJeg - U)c]0 = Up - Uc] 2G = ]2g = , 

(3) 
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where deg is the dipole matrix element. The suscepti- 
bility X can be obtained by considering the steady state 
solution of Eq.ljSJ to the first order of g and write the 
solution as 



P = P" + 
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(4) 



The eg element of p"*" will yield the linear susceptibility x 
at the frequency Up as can be seen by combining Eqs.|(21) 
and l@J: 
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where n is the density of the atoms of the medium. The 
saturation parameter S = 4|GpTiT2 is defined as the ra- 
tio of the control field intensity and the saturation inten- 
sity. The average response of the susceptibility is given 
by 



2Vln2 
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dA, (6) 



where ZJeg is the central frequency of the atomic tran- 
sition |e) < — > \g). Here we consider the frequency of 
the control field ujc is tuned to the line center Ugg- We 
present the behavior of real and imaginary parts of the 
susceptibility as a function of the detuning of the probe 
field in Fig. Q . The real part of susceptibility gives nor- 
mal dispersion. It is clear from Fig. (01) that the slope of 
normal dispersion attains maximum when S '--^ 1 which 



leads to ultra slow light. The imaginary part of (x) ex- 
hibits the absorption dip which becomes deeper with the 
increase in the intensity of the control field as shown in 
Fig. lOs). The spectral width of absorption dip depends 
on the intensity of the control field. This dip is associated 
with coherent population oscillation !l7|. 

In order to compare with experimental results of Baldit 
et al. we need to know the group index Ug which is 
defined by 
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Ug = 1 + 27rWp^^Re(x) 
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where 



(7) 



J 



D 



i{A + i) [52 -I- 2{Sf + + S- A)2(l + iA) + S{i + S- A){-i + 5 + 2f{i + S ~ iS"^ - A) + A)] 
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We denote the integration with respect to A has been 
denoted by (). The unsaturated inhomogeneous absorp- 
tion coefficient of the two level atomic system is defined 
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In the limit of very small detuning of the probe from the 
pump, the analytical expression for the group index for 
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FIG. 2: (a) and (b) The real and imaginary parts respeci 
presence of co-propagating control field G. The common pai 
inhomogeneous absorption coefficient Qfi„h = 6.5 cm~^; inhor 
Ti= 8 ms; transverse relaxation time T2=3 /is. 



inhomogeneous case can be expressed as 



n„ ^ cainhTi 



5(4 + S) 
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It is clear from the above expression that the group in- 
dex varies as 

5-1/2 i^j,gg ^^Yue of S. The group indeji 
attains the maximum value at S* = 0.9. In case of homo- 
geneously broadened two level system the group index is 
given by [2^ 



S 
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where ah — 4:Trujpn\deg\'^T2/ch is the homogeneous ab- 
sorption coefficient. For homogeneous two level system 
the group index varies as at large S and peaks at 
S=0.5. At large S, the group index for a two level sys- 
tem falls much slowly for an inhomogeneous medium as 
compared to the homogeneous case as shown in Fig. ^ . 
We thus find an important difference between inhomoge- 
neously and homogeneously broadened two level systems 
Note that the ratio between inhomogeneous and homo- 
geneous unsaturated absorption coefficient is amh/cth ~ 
^inhT2- The behavior so obtained is consistent with the 
experimental observation. Figure Q shows the varia- 
tion of group index as function of the intensity of the 
control field at different probe detuning. As the detun- 
ing of the probe field is increased the peak of the group 
index shifts toward higher S. The maxima of the group 
index Ug = .65 x 10* for S — 10 Hz occurs at 5 = 1.16 
which corresponds to the group velocity Vg — c/ug — 4.61 
m/s which is higher than what is reported. Further the 



0.08 



0.06 -. 



0.04 




0.02 - 



FIG. 3; Variation of the term in squared bracket of Eqs. HlOjl 
and mil as a function of intensity of the control field for 
inhomogeneous and homogeneous cases of two level system. 



to understand the nature of the experimental results of 
Baldit et al. However the results derived here are appli- 
cable to any system which can be modelled by inhomo- 
geneously broadened two level system. 
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FIG. 5; Group index variation with the detuning of the probe 
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